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Abstract 

It is shown that the Dirac sea can be uniquely defined for the Dirac 
equation with general interaction, if we impose a causality condition 
on the Dirac sea. We derive an explicit formula for the Dirac sea in 
terms of a power series in the bosonic potentials. 

The construction is extended to systems of Dirac seas. If the system 
contains chiral fermions, the causality condition yields a restriction for 
the bosonic potentials. 



1 Introduction 



The Dirac equation has solutions of negative energy, which have no mean- 
ingful physical interpretation. This popular problem of relativistic quantum 
mechanics was originally solved by Dirac's concept that all negative-energy 
states are occupied in the vacuum forming the so-called Dirac sea. Fermions 
and anti-fermions are then described by positive-energy states and "holes" in 
the Dirac sea, respectively. Although this vivid picture of a sea of interacting 
particles is nowadays often considered not to be taken too literally, the con- 
struction of the Dirac sea also plays a crucial role in quantum field theory. 
There it corresponds to the formal exchanging of creation and annihilation 
operators for the negative-energy states of the free field theory. 



e-print archive: |http://xxx. lanl.gov/abs/hep-th/9765006] 
^Supported by the Deutsche Forschungsgemeinschaft, Bonn. 



964 



DEFINITION OF THE DIRAC SEA 



Usually, the Dirac sea is only constructed in the vacuum. This is often 
considered to be sufficient, because the interacting system can be described 
by a perturbation of the vacuum. 

Unfortunately, the situation is more difficult: In relativistic quantum 
mechanics with interaction, the fermionic wave functions are solutions of 
the Dirac equation 

{i^ + B-m)^ = 0, (1.1) 

where the operator B is composed of the bosonic potentials (for example, 
we can describe the electromagnetic interaction by choosing B = eJ^ with 
the electromagnetic potential A). 

In contrast to the free Dirac equation (i^ — m) ^' = 0, it is not obvious 
how to characterize the negative-energy solutions of (11. ip . Qualitatively, the 
problem is that the perturbation B leads to a mixing of the free solutions and 
destroys the natural splitting into solutions of positive and negative energy. 
As a consequence, it is not clear how the Dirac sea of the system (II. ID can 
be constructed. We point out that this problem is not solved by a simple 
perturbation expansion in B; it is then hidden in the non-uniqueness of this 
expansion (see section [2] for details). In quantum field theory, the problem 
of defining the Dirac sea is even more complicated, because the virtual pair 
creation/annihilation must be taken into account. We will not deal these 
problems here and restrict to the limit of "classical" potentials and wave 
functions. Nevertheless, our considerations are also relevant for quantum 
field theory, because it is in many situations (e.g. for a quantum system in 
a classical background field) preferable to use the Dirac equation (jl.ip as 
the starting point for the fermionic field quantization. In this sense, the 
construction of the Dirac sea of ()l.ip is preliminary for the description of 
interacting quantum fields. 

We conclude that the definition of the Dirac sea is basic for a reasonable 
physical interpretation of the Dirac equation (jl.ip . In the present paper, we 
will discuss the difficulty in constructing the Dirac sea and finally solve the 
problem in terms of a formal perturbation expansion in B. Before starting 
the analysis, we describe the problem in more mathematical terms: Every 
solution of the free Dirac equation (z^ — m) ^' = is a linear combination of 
plane wave solutions of the form 

with a 4-spinor % which is independent of t and x. The sign of lo gives 
a natural splitting of the solutions into solutions of positive and negative 
frequency. Identifying frequency and energy via Planck's formula, these so- 
lutions are commonly called the positive and negative energy solutions of the 
free Dirac equation. Since the simple identification of frequency and energy 
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might lead to confusion (sometimes the "energy" of a negative-frequency 
state denotes the positive energy of the corresponding anti-particle state), 
we prefer the notion of positive and negative "frequency" in the following. 
We denote the negative-frequency solutions by where k is the momen- 
tum and a = 1, 2 are the two spin states (for an explicit formula for ^'^^ see 
e.g. [U). If the states ^'^^ were normalized with respect to the usual scalar 
product 



(^1$) = / {m ^){t,x) dx , ^ = ^*70 , (1.2) 

we could form the projector P<.i,^ > on the one-dimensional subspace <^'r: > 

ka rCft 

by 

(^<*.^.> ^) ^) = (^fca(*' ^) ^^t^) 7° ^(i, y) dy 

In this sense, the product ^'g^(x) ^'^^(y) would be the kernel of the projector 
on <\I'^^>, and the sum over all negative- frequency states would yield the 
projector on the whole Dirac sea. Unfortunately, the wave functions 
are not normalizable. We could arrange normalizable states by considering 
the system in finite three-volume, but we do not want to do this here. It 
is more appropriate for our purpose to formally build up a projector on all 
negative-frequency states by integrating over the momentum parameter 

P{x,y) = Y: L ' (1-3) 

a=l,2 

which can be rewritten as the integral over the lower mass shell 



/ir3 (27r) 

(G denotes the Heavyside function Q{x) = 1 for x > and Q{x) = other- 
wise). P{x,y) is a well-defined tempered distribution which solves the free 
Dirac equation {i^^ — m)P{x,y) = 0. We can use it to characterize the Dirac 
sea in the vacuum. Our aim is to introduce a corresponding distribution P 
for the Dirac equation with interaction (jl.ip . The construction of P must 
be unique in a sense which we will discuss and specify later. We will as- 
sume the perturbation ;B to be a differential operator on the wave functions. 
Furthermore, it shall be Hermitian with respect to the (indefinite) scalar 
product 

<^ [ = j ^{x) (^{x) d'^x . (1.4) 

For an electromagnetic potential B = e^, these assumptions are satisfied 
because 4 = ^^4^^^- In addition, B can be composed of the scalar, pseu- 
doscalar, pseudovector and bilinear potentials as e.g. discussed in [2j. Ac- 
cording to [3], i3 also allows for the description of the gravitational field. 
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2 Non-Uniqueness of the Simple Perturbation Ex- 
pansion 

Our first idea for the construction of P is to calculate solutions of (11. ip 
with a perturbation expansion in B and to define P in analogy to (|1.3p by 



P{x,y) = E L^kai^)^kaiy)dk . (2.1) 

We start with a discussion of this method in a perturbation calculation 
to first order. This is quite elementary and will nevertheless explain the 
basic difficulty. For the perturbation calculation, we need a Green's function 
s{x, y) of the free Dirac operator, which is characterized by the distributional 
equation 

{i^.^-m)s{x,y) = 5\x-y) . (2.2) 

To first order, the perturbed eigenstates ^'g^ are then given by 

^kS^) = *fca(^) - / dSs{x,y)By^^^^{y) + 0{B^) , (2.3) 

as can be verified by substituting into (jl.ip . We insert this formula into 
(12. ip and obtain 

P{x,y) = P{x,y)-J d^z [s{x,z)B,P{z,y) + P{x,z) B, s*{z,y)]+0{B^), 

(2.4) 

where we used that B is Hermitian with respect to the scalar product (jl.4p . 
and where s*{z,y) is given by s*{z,y) = s(7/, z)^ 7^^. It is convenient to 
view the distributions s{x,y), P{x,y) as integral kernels of corresponding 
operators s,P. Then we can write (|2.4p with operator products 



P = P - sBP - PBs* + OiB"^) , (2.5) 

where the superscript '*' denotes the adjoint with respect to the scalar prod- 
uct (nr 



Equation (j2.5p gives a possible definition for P. As apparent problem, the 
construction depends on the choice of the Green's function. For example, we 
could have chosen for s either the advanced or the retarded Green's function 
^rn^^mi which are in momentum space as usual given by 

vri\ r ^ + m 1^ + m 



(2.6) 

More systematically, the arbitrariness of our construction is described as 
follows: According to p.2p . the difference between two Green's functions is 
a solution of the free Dirac equation. We can thus represent s in the form 

s{x,y) = s^ni{x,y) + a{x,y) 
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where a[x^ y) is in the x-variable a hnear combination of the plane-wave 
solutions, i.e. 

a{x,y) = L^kai^)c^,{y) dk 

with (generally complex) functions c^^{y), where = 3,4 denote the 

plane-wave solutions of positive frequency. We substitute into (|2.5p and 
obtain 

P = p - sl^BP - PBs^, - {aBP + PBa*) + OiB"^) . (2.7) 



The expression in the brackets maps solutions of the free Dirac equation into 
each other and vanishes otherwise. We can thus write it in the form 




(2.8) 

with suitable functions gab{ki, ^2)- This representation of P can also be un- 
derstood directly: The contribution (j2.8p describes a mixing of the solutions 
^'^^ of the free Dirac equation. To the considered first order in B, it vanishes 
in the Dirac equation (i^ + B — m) P = 0. Thus we cannot fix this contri- 
bution with the Dirac equation, it remains undetermined in our method. 
According to (12. 7p . this is the only arbitrariness of the construction; the 
other contributions to P are unique. 

In higher order perturbation theory, the non-uniqueness can be under- 
stood similarly, although the situation is more complicated: For a given 
Green's function s, we can construct a solution of the Dirac equation 
(jl.ip by the formal perturbation series 

00 

= E(--^r^.^a > (2-9) 

n=0 

as is verified by substituting into (jl.ip . Actually, this is a very special 
ansatz. For example, we can use different Green's functions in every order 
of the perturbation calculation, which leads to the more general formula 

00 

^ka = ^ka + ^^"^ • • • ^^^^ ^^^^ *fea (^-^O) 

n=l 

with a whole series of arbitrary Green's functions s^^\ s^^\ etc.. Once we 
have a formula for the non- uniqueness of P can again be discussed by 
substituting into (|2.ip . In generalization of ()2.8[) . the arbitrariness of the 
construction is described by a contribution to P{x, y) of the form 
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which mixes perturbed eigenstates ^'g^ and vanishes in the Dirac equation 
+ B — m) P = 0. The dependence of Qabiki, on B and on the Green's 
functions s^"^ is rather involved, however, and we need not go into the details 
here. 

To summarize, a simple perturbation expansion in B is not unique and 
therefore does not allow a meaningful definition of P. In the ansatz (I2.10p . 
for example, we should find a way to specify the Green's functions s^"^. This 
cannot be done with the Dirac equation (II. ip . and we must therefore look 
for additional input to completely determine P. Our basic idea is to apply 
some causality principle. For example, it might seem a reasonable condition 
to impose that P{x,y) only depends on B in the "diamond" {L]^ n Ly) U 
(L^nLp, where 

= {y| >0, yO-xO >0} (2.11) 

= {y| >0, y°-xO <0} (2.12) 

denote the future and past light cones around x, respectively. If we want 
to study conditions of this type, it is no longer useful to look at the pertur- 
bation expansion for the individual states ^g^(x) (because these states only 
depend on one argument x). We must take into account for the perturbation 
expansion that P is composed of many states in a specific way. 



3 The Causal Perturbation Expansion 

In preparation, we first describe how the perturbation expansion for the 
advanced and retarded Green's functions can be performed uniquely: The 
support of the distribution s^{x,y) is in the future light cone y € (this 
can be checked by calculating the Fourier transform of ()2.6p with contour 
integrals). As a consequence, the perturbation operator B{z) only enters 
into the operator product 

{staBsl^){x,y) = d^z sl^{x,z) B{z) s]i,{z,y) (3.1) 



for z e L^nL^. In this sense, the expression (13. ip is causal. Especially, the 
support of (|3.ip is again in the future light cone. It follows by iteration that 
the higher powers 

S-m ^ SjYi B • • • B B 

are also causal and have their support in the upper light cone. We define the 
perturbed advanced Green's function as the formal sum over these operator 
products, 

oo 

sl = Y.i-'l^^f ■ (3-2) 

k=0 
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Accordingly, the perturbed retarded Green's function is defined by 

oo ^ 

— ^ {~^m, ^) ■ (3-3) 
k=Q 

These operators satisfy the defining equations for the perturbed Green's 
functions 

{i^-m + B)s'^ = 1 = {i^-m + B)s^ , (3.4) 

as is verified directly. 

Notice that the perturbation expansion for the Green's functions becomes 
unique by the condition that the contribution to s^, to every order has 
its support in the future and past light cones, respectively. We want to 
take this construction as the guiding line for the perturbation expansion of 
P. Unfortunately, the method cannot be directly applied to the Dirac sea, 
because the distribution P{x, y) does not vanish for space-like y — x, and we 
thus have no notion of causality. As way out, we decompose the free Dirac 
sea in the form ^ 

P{x,y) = ^{Pm{x,y) - km{x,y)) (3.5) 
with the tempered distributions 

Pm{x,y) = J ^{f^ + m)6{k^-m^)e'''''^--y^ (3.6) 

^ (^ + m) 5{e - m") e(fcO) e^^^^-^^) (3.7) 

(e denotes the step function e{x) = 1 for x > and e{x) = — 1 otherwise). 
We also consider pmix,y) and km{x,y) as integral kernels of corresponding 
operators pm, km- The operator pm is built up as a formal sum over the 
projectors on all solutions of the Dirac equation and can be viewed as a 
spectral projector of the free Dirac operator. The definition of km differs 
from Pm by a relative minus sign for the states on the upper and lower mass 
shell. As a consequence of this relative minus sign, the Fourier integral (13. 7p 
vanishes if y — x is space-like (this can be seen from Lorentzian invariance 
and a symmetry argument for k = {0,k)). Thus km{x,y) is causal in the 
sense that it has the support in the light cone y € U L^. This makes 
it possible to uniquely express its perturbation expansion in terms of the 
perturbed Green's functions: We substitute the distributional equation 

lim ( — ^ = 27rz 6(x) 

o<e-^o \x -ie x + iej 

into the formula for km in momentum space, 
km{p) = {'^ + m) 5{p^ - w?) e{p^) 
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id + m) lim 

("d + m) lim 



p2 _ _ ^£ 



p2 _ _ jgpO 



— + ie 
1 



p2 _ j^2 _|_ j^^pO 



and obtain with (j2.6p a simple relation between fc^ and s^, s^, 



2TTi ^ ™ 



(3.8) 



We extend this relation to the case with external fields: 



Def. 3.1 We define the operator km by 
with the Green's functions 13.^) . [373\) . 

According to (j3.4p . fc^ really is a solution of the Dirac equation (i^ + B — 
m) km = 0. 

In order to explain the significance of this construction, we point out 
that the factor e{k^) in (j3.7|) describes the splitting of the solutions of the 
free Dirac equation into solutions of positive and negative frequency. With 
the introduction of km, we were able to uniquely generalize this splitting to 
the case with external fields. This solves the basic problem in defining the 
Dirac sea. It remains to perform the perturbation expansion for pm- On a 
formal level, this is very easy, because we can remove the relative minus sign 
for the positive and negative frequency states by taking the absolute value 
of km, 

~ formally [^7 in\ 

Vm := • (3-10) 

This gives a unique definition for pm- Since km is composed of eigenstates 
of the perturbed Dirac operator with eigenvalue m, it follows automatically 
that (i(^ + B — m) pm = 0. 

Unfortunately, it requires some effort to convert the formal relation (j3.10p 
into a mathematical definition. The problem is that the square km is ill- 
defined; furthermore we want to write pm as a power series in B. These 
problems are solved in the following theorem. The reader who is not so in- 
terested in the technical details and the combinatorics of the expansion may 
skip the proof. For the statement of the theorem, we need some notation: 
We work with the Green's function 

Sm = lis^m + S^k) , (3.11) 
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which has the advantage of being Hermitian (with respect to the scalar 
product (II. 4p ). Furthermore, we introduce the series of operator products 

oo oo oo 

k=0 k=0 k=0 



and set for Q C N 

Fm{Q,n) 



Pm iin e Q 
km if n ^ Q 



Theorem 3.2 The relations \3. 9^) . [3A^) uniquely determine the perturba- 
tion expansion for km and pm ■ We have the explicit formulas 

oo 

~km = Y.^-iT^f^ h<km{bmkm?^b'^ (3.12) 
0=0 

oo L 2 J 
/3=0 a=0 

with the coefficients 

c(0,0) = 1 , (3.14) 



c(a,/3) 



•—a+l \ / 



and the operator products 

Gm{0C^f3) = E {-iTTf^b<FmiQ,l)bmkmbmFm{Q,2) 

QeV{l3+l), #Q=2a+l 

X bmkmbm ' ' ' bmkmbm FmiQj + 1) j (3.16) 

where V{n) denotes the set of subsets of {1, . . . ,n} (we use the convention 
ni = 1 forl< 0). 

Proof: Notice that (i^ + B — m) bm = 0. Since all operator products in 
(j3.12p . (j3.16p have a factor bm at the left, the operators pm, km are solutions 
of the Dirac equation 

{i^ + B — m) Pm = = {i^ + B — m) km 

Thus the theorem gives a possible perturbation expansion for pm and km- We 
must verify that the conditions ()3.9p . (l3.10p are satisfied and show uniqueness. 



According to (13. Sh . the advanced and retarded Green's function can be 
written in the form 

~ ^'^ km j Sm ~ ^rn ~ i^^ km ■ (^•1''') 
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We substitute the sums ([32]), ([23]) into dSlQ]) . 

1 oo 

^™ = ^ E ((-«^ «™ - (-^- ' (3-18) 

insert (13.171) and expand. This gives a sum of operator products of the form 

CiB C2B ■ ■ ■ B with Cj = km or Cj = Sm 

The contributions with an even number of factors km have the same sign 
for the advanced and retarded Green's function and cancel in ()3.18p . The 
contributions with an odd number of fem's occur in every Green's function 
exactly once and have opposite sign. Using the notation 

we can thus rewrite (|3.18p in the form 

00 

km = E(-l)' E (-)*^-' 

1=0 QeV(l+l), #Qodd 

X CmiQ, 1) B CmiQ, 2)B---B CmiQ, I) B CmiQ, I + 1) 

After reordering the sums, this coincides with (j3.12p . 

Next we want to give the relation (j3.10p a mathematical sense. For this, 
we consider m > as a variable mass parameter. Then we can form products 
of the operators pm, km by manipulating the arguments of the distributions 
in momentum space. For example, we have with (j3.6p 

Pm{k)Pm'{k) = + m) 6{k^ - m^) + m') 5{k^ - (m'f) 

= {k^ + {m + m')^ + mm!) 5{m? — (m')^) 5{k'^ — m?) 

= {k^ + (m + m')^ + mm') - — 6{m — m') 6{k^ — m?) 

= 5{m-m')pm{k) , (3.19) 

and similarly with (j3.7l) . 

Pmkm' = km' Pm = 5{m - m') km (3.20) 
kmkm' = 5{m-m')pm ■ (3-21) 

This formalism has some similarity with the bra/ket notation in quantum 
mechanics, if the position variable x is replaced by the mass parameter m. 
Equation (|3.19p can be understood directly from the fact that pm are the 
spectral projectors of the free Dirac operator; the relations (j3.20p . (j3.2ip 
reflect the relative minus sign in km for the states on the upper and lower 
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mass shell. Especially 

= Pm Pm'- This relation can be 

extended to the case with interaction, 

Pm Pm' ~ f^m ^m' : (3.22) 

and gives a meaningful square of (I3.10p (we will see in a moment that km km' 
vanishes for m ^ m'). If our construction ensures that pm is a positive 
operator, (j3.22p is even equivalent to (j3.10p . 

We calculate the product km km' explicitly: The definitions (j3.6p . (j3.7p 
and (jXTTIl . ffZB yield in analogy to (f3l9]) the formulati 



Pm Sm' — Sm' Pm — PP ( , I Pm (3.23) 

\m — m' / 

^m ^m' ^m' ^m PP I - | km (3.24) 

ym — m' / 

Sm Sm' — PP ( 7 I (Sm ~ Sm') i (3.25) 

V m — m' / 

where PP(x~-'^) = ^ limo<e_to[(2^ + i^)^^ + {x — ie)'^] denotes the principal 
value. As a consequence, the operator products with factors Sm,Sm' are 
telescopic, i.e. 

n 

km {B SmT {sm' S)""^ ^m' =0 forn > 1. (3.26) 

This allows us to easily carry out the product bm in the expression 

km b^' km' = S{m - m!) pm ■ (3.27) 

With this formula, we can calculate the square of ()3.12p to 

oo 

km km' = 5{m - m') ^ {-iT,fP^+^P^ b< {km bmf^' Pm {bm kmf^' b> . 
I3i,f32=0 

(3.28) 

We could continue the proof by verifying explicitly that the product 
Pm Pm' with Pm according to (I3.13P coincides with (j3.28p . This is a straight- 
forward computation, but it is rather lengthy and not very instructive. We 
prefer to describe how the operator products (j3.16p and the coefficients p.lSp 
can be derived. In order to keep the proof better readable, we make some 
simplifications: Since the factors b^, b^ cancel similar to (j3.27p in telescopic 
sums, we can omit them in all formulas without changing the multiplica- 
tion rules for the operator products. Then all operator products have km or 
Pm as their first and last factor, and we can multiply them with the rules 



^Note added in January 2009: As noticed by A. Grotz, in (I3.25P the summand 7r'^(5(m — 
m') Pm is missing. This error is corrected in the recent paper arXiv:0901.0334 [math-ph]. 
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(j3.19p . (j3.20p . and (j3.2ip . Since all these rules give a factor 5{m — m'), we 
will in any case get the prefactor 6{m — m') in (I3.28p . Therefore we can 
just forget about all factors 6{m — m') and consider all expressions at the 
same value of m. Furthermore, we will omit the subscript 'm' and write the 
intermediate factors 6 as a dot After these simplifications, we end up 
with formal products of the form 

Fi . F2 . F3 . • • • .Fn with Fj = k or Fj =p (3.29) 

and have the multiplication rules 

p"^ = k'^ = I ^ pk = kp = k . (3.30) 

We must find a positive operator p being a formal sum over operator products 
([3:29]) such that 

00 

p^ = {-i7,f''-+^P'^{k.f^^p{.kff'^ . (3.31) 

/3i,/32=0 

In this way, we have reduced our problem to the combinatorics of the opera- 
tor products. As soon as we have found a solution p of (j3.3ip . the expression 
for pm is obtained by adding the subscripts 'm' and by inserting the factors 
^m, Relation (j3.22p follows as an immediate consequence of (|3.3ip . 

The basic step for the calculation of p is to rewrite (I3.3ip in the form 
p2 = p+A with A = {-i^?^'^^^'{k.f^'p{.kf^^ . 

(/3i,/32)^(0,0) 

(3.32) 

The operator p is idempotent and acts as the identity on A, Ap = pA = 
A. Therefore we can take the square root oi p + A with a formal Taylor 
expansion, 

P = = P + £ (-1)'^"'' ^^^r^ A- , (3.33) 

n=l 

which uniquely defines p as a positive operator. 

It remains to calculate A". If we take the nth power of the sum in 
(j3.32p and expand, we end up with one sum over more complicated operator 
products. We first consider how these operator products look like: The 
operator products in (j3.32p all contain an even number of factors k and 
exactly one factor p. The factor p can be the 1st, 3rd,. . . factor of the 
product. Each combination of this type occurs in A exactly once. If we 
multiply n such terms, the resulting operator product consists of a total odd 
number of factors p,k. It may contain several factors p, which all occur at 
odd positions in the product. Furthermore, the total number of factors p 
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is odd, as one sees inductively. We conclude that consists of a sum of 
operator products of the form 

{k . k .y^ p.k.{k.k .y^ p.k.{k.k .y^ ■■■ {k.k .)^'"+' p{.k. ky^-+^ (3.34) 

with a, qj > 0. We set /3 = 2a + J2j Qj- Notice that the number of factors p 
in ()3.34p is 2q + 1; the total number of factors p,k is 2/? + 1. The form of 
the operator product gives the only restriction < 2a < /3 for the choice of 
the parameters a, /3. 

Next we count how often every operator product ()3.34p occurs in the 
sum: The easiest way to realize (j3.34p is to form the product of the a + 1 
factors 

'{k.k.y^ p {.k.ky^+^] Uk.k.y^^^ p {.k.ky^^^] ■ ■ ■ \{k.k.y^"+'+^ p {.k.kY^--+^ 

(3.35) 

However, this is not the only possibility to factorize (j3.34p . More precisely, 
we can apply to each factor in (j3.35p the identities 

{k.k.yp{.k. kf = [{k.k .y p] [p{.k. ky] 
{k.k.yp{.k.kY = [{k.k.fp] [{k .k .y-' p{.k .kf] 
{k .k .y p{.k .kf = [{k .k .y p{.k .kf-'] [p{.k.kY] 

By iteratively substituting these identities into ()3.35p . we can realize every 
factorization of ()3.34p . Each substitution step increases the number of factors 
by one. The steps are independent in the sense that we can fix at the 
beginning at which positions in (j3.35p the product shall be split up, and 
can then apply the steps in arbitrary order. There are (a + 1) + {qi — 1) + 
Z]j=2^"^ + {Q.2a+2 — 1) = /? — (a + 1) positions in (I3.35P where we could split 
up the product (in the case gi = or q2a+2 = 0, the counting of the positions 
is slightly different, but yields the same result). Since we want to have n 
factors at the end, we must choose n — (a+1) of these positions, which is only 
possible for a + 1 < n < /? and then gives (/? — q — l)!/((n — a — 1)! (/? — n)!) 
possibilities. 

Combining these combinatorial factors with the constraints < 2a < /?, 
a + 1 < n < fi gives for n > 1 

oo min(n-l,[f]) . V 
j3=n a=0 \ I Qgp(/3+l), #Q=2a+l 

X {-iTrf^ F{Q, I) .k. F{Q, 2) . k . ■■■ .k. F{Q, 13 + 1) (3.36) 
with -F((5, n) = p for n ^ Q and -F(Q, n) = k otherwise. Notice that the last 



sum in (|3.36p runs over all possible configurations of the factors p, k in the 
operator product (j3.34p for fixed a, [3. We finally substitute this formula into 
(I3.33P and pull the sums over a, (3 outside. This gives the desired formula 
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for p. ■ 

In order to illustrate the derived formulas for p and A;, we give the contribu- 
tion up to third order in more detail: 

kxYi — kffi kffi 13 Sffi SfYi 13 kfn 

13 kffi 13 Sfn 13 Sfn km 13 Sm 13 Sm 13 Sm 
~\-7r Sm 13 km ^ kfn ^ km ~l~ km ^ ^ ^ km 

+7r^ km B km B Sm B km + tt"^ km, B km B km, B Sm + 0{B'^) 

Pm — Pm Pm B Sm Sm B Pm 

~^^m B Sm B Pm Sm B Pm B Sm ~1~ Pm 

B SmB 

^ Pm B km B km km B km B Pm 

B SmB Sm.B 

Pm ^m 

B SmB 

Pm ^ ^m 

i~' Pm 

B SmB 

^m Pm B Sm B Sm B Sm 
+ — (s^ B Pm, B kjn B kfn + Pm B Sm, B km B km, 



~^Pm B km B Sm, B km, + Pm B km B km B Sm) 



~\~km 13 kfn 13 Sm 13 Pm ~\~ km ^ km 13 Pm 

B Sm) + OiB^) 

The theorem gives precise formulas for the perturbation expansion of the 
Dirac sea. Both the combinatorics of the factors km,Pm, Sm and the numeri- 
cal prefactors are a non-trivial result and, as far as the author knows, cannot 
be understood intuitively. 

We call the perturbation expansion of this theorem the causal perturba- 
tion expansion. It allows to uniquely define the Dirac sea by 

P{x,y) = ^{Pm-km){x,y) 



4 Generalization to Systems of Dirac Seas 

In the previous section, we defined the Dirac sea for a system of interacting 
fermions of mass m. A realistic model, however, is composed of several types 
of fermionic particles with masses mi, . . . ,mf. Furthermore, the fermions 
of zero mass may (like the neutrinos in the standard model) occur only as 
left or right handed particles. The perturbation B will in general mix up the 
eigenstates to different masses and will in this way describe an interaction of 
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all the fermions. We will now extend the previous construction to this more 
general setting. 

First we must generalize (|1.3p and define a distribution y) which 
describes the system in the vacuum: In order to distinguish the chirality 
of the zero-mass fermions, we introduce (4 x 4)-matrices Xi, . . . ,X/. For 
the zero-mass fermions nij = 0, they can be either Xj = 1, Xj = xl, or 
Xj = xr, where Xl/r = ^(1^7^) are the chiral projectors. For rrij ^ 0, they 
must coincide with the identity Xj = 1. The Dirac seas of the individual 
types of fermions are then described by Xj ^ {pm^ — km.j ){x,y)- The remaining 
question is how to build up P{x,y) from the individual Dirac seas. In view 
of the configuration and the interactions of the fermions in the standard 
model, one might want to use combinations of sums and direct sums 

(e.g. with ? = 1, . . . , 8 = 2(3 + 1) running over the color, lepton, and isospin 
index, and with the index a = 1,...,3 to distinguish the three fermion 
families. It seems reasonable to use the ordinary sum over a because the 
families show the same interactions). From the mathematical point of view, 
however, it is easier to use only direct sums 

^ 1 

P{x,y) = X; -(p„, - k^,){x,y) . (4.2) 

1=1 ^ 

This is no loss of generality, because the ansatz (14. Ih can be obtained from 
()4.2p by taking a suitable partial trace over the /-index (in our example, 
by choosing / = 24 = 3 • 8 and forming the trace over the three families). 
For the perturbation expansion, we can also restrict ourselves to the ansatz 
(j4.2p . because the perturbation expansion for ()4.ip is obtained by taking the 
partial trace of P(x, y) (see [4j for a more detailed discussion of this method). 
Therefore we must in the following only consider a P{x, y) of the form (j4.2p : 
it is called the fermionic projector of the vacuum. 

It is convenient to use a matrix notation in the direct sum: We set 

p{x,y) = ®i^iPmi{x,y) , 

and define the matrices 

/ 

X = , Y 

1=1 

which are called chiral asymmetry matrix and mass matrix, respectively (m 
is an arbitrary mass parameter; e.g. one can choose m = max^mj). Then 
we can write the fermionic projector as 

P{x,y) = X^{p{x,y)-k{x,y)) . (4.3) 



k{x,y) = ef^^ km,{x,y) 

1 ^ 
= — 
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Since rrij = for X 7^ 1 and because Pm=o, km=o anti-commute with 7^, we 
have alternatively 

P{x,y) = ^{p{x,y)-k{x,y))X* , (4.4) 

where X* = 'y^X^^^ is the adjoint with respect to the scalar product (jl.4p . 
The fermionic projector is a solution of the free Dirac equation 

{i^^-mY) P{x,y) = 

In order to describe the interacting system, we again insert a differential 
operator B into the Dirac equation. Thus the fermionic projector P{x, y) is 
supposed to be a solution of the Dirac equation 

(i^^ + B-mY) P{x,y) = . (4.5) 

B may be non-diagonal in the "Dirac sea index" /; we assume it to be Her- 
mitian with respect to the scalar product 

/ r 

<^ I = ^ / ^'^(x) $/(x) d'^x 
1=1 



The perturbation expansion for k and p can be carried out exactly as in 
the previous section: We define the advanced and retarded Green's functions 
by 

s^{x,y) = ej=i 5^^.(2;, y) , s^{x,y) = ej=i s^^.(x,y) 

Their perturbation expansion is, in analogy to (j3.2p . (|3.3p . uniquely given by 

00 00 
gv ^ Y.i^-s^Bfs' , 5^ = Y,{-s^Bfs'^ . (4.6) 

A:=0 k=0 



Theorem 4.1 The perturbation expansion for p and k is uniquely deter- 
mined by the conditions 

1 ~A\ ~ formally 



k = -^,{^-~s^) , p^^"^-"y^k2 , (4,7) 



We have the explicit formulas 



00 L 2 J 



~k = Y.(-iT:fl^ b< k{hkfH> , p = ^ ^c(a,/?)G(Q,/3) 

/3=0 /3=0o=0 
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with 

c(0,0) = 1 

t a^ 4^ t (2n-3)!! ( (3 - a - l\ r a ^ a 

c(a,B) = > (-1)"+^ ^ ; — I forB>l and 

n=a+l \ / 

Gif,g) = J2 {-iTi?^h< F{Q,l)bkhF{Q,2)bkh ■■■ 

Q&V{I3+1), #Q=2a+l 



xbkbF{Q,p + l)b> , 
where V{n) is the set of subsets o/{l, . . . ,n} and where we used the notation 

oo oo oo 

b< = Y.i-'^)' ' b = J2i-Bs)>'B , b> = 

k=0 k=0 k=0 



Proof: Follows exactly as Theorem [X 



After this straightforward generalization, we come to the more interesting 
question of how P can be defined. Our first idea is to set in generalization 
of ^ 

Pix,y) = X^{p-k){x,y) . (4.8) 

This is not convincing, however, because we could just as well have defined 
P{x,y) in analogy to (|4.4|) by P = ^{p — k) X* , which does not coincide 
with (j4.8p as soon as X,X* do not commute with B. It turns out that 
this arbitrariness in defining the Dirac sea reflects a basic problem of the 
causal perturbation expansion for systems with chiral asymmetry. In order to 
describe the problem in more detail, we consider the perturbation calculation 
for k to first order: According to (|4.6p . (j4.7p . we have 

k = k - — {s^ B s"" - s"" B s"") + 0{B'^) (4.9) 
27ri 

= k - sBk - kBs + 0{B^) 

This expansion is causal in the sense that k{x,y) only depends on B in the 
"diamond" (L^ n L^) U (L^ n L^), as is obvious in (fO|) . It is not clear, 
however, how to insert the chiral asymmetry matrix into this formula. It 
seems most natural to replace all factors k by Xk, 

(Xk) = Xk - sBXk - XkBs + 0{B^) . (4.10) 

This formula really gives a possible perturbation expansion for the system of 
Dirac seas. Unfortunately, it cannot be expressed similar to (j4.9p with the 
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advanced and retarded Green's functions, which means that the causaHty of 
the expansion is in general lost. In order to avoid this problem, one might 
want to insert X at every factor s, /c, 

{Xk) = Xk - XsBXk - XkBXs + 0{B^) 

= Xk ^—iXs'BXs'^-Xs'^BXs'^) + 0{B'^) (4.11) 

zvrz 

Similar to (j4.9p . this expansion is causal. In general, however, it does not 
give a solution of the Dirac equation (z^ + B — m) k = which does not 
make sense. 

The only way out of this dilemma is to impose that the perturbation 
expansions ()4.10p and (j4.1ip must coincide. This yields a condition on the 
perturbation operator B, which can be characterized as follows: We demand 
that 

Xs^BXs'^ = BXs^ = Xs^Bs^ . (4.12) 

Since the operator s^^q anti-commutes with 7^, we have Xs^ = X* . 
Substituting into the second equation of (I4.12P yields the condition X* B = 
B X. Since X is idempotent, this condition automatically implies the first 
equation of (j4.12p . We formulate the derived condition for the whole Dirac 
operator i(^ + B — mY and thus combine it with the fact that chiral fermions 
are massless (i.e. X*Y = YX = Y) and that X is composed of chiral 
projectors (which implies that X*^ = (^X). 

Def. 4.2 The Dirac operator is called causality compatible with X if 

X* {i^ + B- mY) = {i^ + B- niY) X . (4.13) 

In the perturbation expansion to higher order, the condition (j4.13p allows 
to commute X through all operator products. Using idempotence = X, 
we can moreover add factors X to the product, especially 

X CiBCiB ■■■ BCn = XCi B XCi B ■■■ B XCn 

with 

Cj = p, Cj = k ov Cj = s . 

This ensures the causality of the perturbation expansion. For a Dirac op- 
erator which is causality compatible with X, the fermionic projector in the 
external field is uniquely defined by ()4.8p . 

5 Discussion, Outlook 

In this paper, we gave the formal definition of the Dirac sea in the presence 
of external fields. The method differs considerably from earlier attempts 
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to solve the external field problem (see e.g. [2] and the references therein). 
Namely, in these previous approaches, the Dirac sea was always constructed 
as the "negative frequency solutions" of the Dirac equation. The basic prob- 
lem of this concept is that the notions of "positive" and "negative" frequency 
do not make sense in the case with general interaction. Therefore the con- 
struction was always limited to potentials which are either static or have 
an only adiabatic time dependence. As shown in this paper, the notion of 
"negative frequency states" is not essential for the description of the Dirac 
sea. For a general definition of the Dirac sea, it must be given up and must 
be replaced by a suitable causality condition. In this way, it becomes possi- 
ble to define the Dirac sea in the presence of potentials with arbitrary time 
dependence. Although the details of the perturbation expansion are a bit 
complicated, the basic concept is very simple. The construction is explicitly 
covariant. It puts the usual "hole" -interpretation of the Dirac equation on 
a satisfying theoretical basis. 

In order to clarify the connection to the usual definition of the Dirac 
sea, we describe how our definition simplifies in the limit of static potentials: 
If considered as multiplication operators, static potentials map functions of 
positive (negative) frequency into functions of positive (negative) frequency. 
Since p, k, and s are diagonal in momentum space, they clearly also preserve 
the sign of the frequency. Thus we have 



quency, respectively. The operators p and k only differ by a relative minus 
sign for the states of positive and negative frequency, 



Using this relation together with (15. 1|) . we can replace pairs of factors p by 
pairs of factors k, e.g. 



■■■pB---pB--- = ■■■pB---pB--- {P+ + P-) 

= P+(- --kB-'-kB---) + p-{--- i-k) B--- i-k) B---) 
= ■■■kB---kB--- 



where the dots '• • •' denote any combination of the operators s, k, p, and 
B. This allows us to simplify the formula for p by only using exactly one 
factor p in every operator product. After going through the details of the 
combinatorics, one obtains the formula 





±P^k 



oo 



p = Y,{-i7:fH< p{hkfH 



6=0 

Thus the Dirac sea (14. 8 p can be written as 



P{x,y) = Y.{-infH< -X{p-k) {bkf'b 
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This equation shows that P(x, y) is composed of the negative-frequency 
eigenstates of the Dirac operator (notice that the expression in the brackets 
[• • •] is the fermionic projector of the vacuum and that ah other factors pre- 
serve the sign of the frequency). Thus, for static potentials, our definition is 
equivalent to the usual concept of "negative frequency states." On the other 
hand, this consideration illustrates in which way our definition goes beyond 
the usual picture. 

In order to get a better understanding of the time-dependent situation, 
we next consider a scattering process. For simplicity, we use the elementary 
framework of [IJ, but our consideration also applies to the operator algebra 
and Fock space formalism as e.g. described in We first recall how a 
scattering process is commonly described in the classical Dirac theory. We 
assume the scattering to take place in finite time to < ^ < ^i- This means 
that the external perturbation B in (jl.ip vanishes outside this time interval, 

B{t,x) =0 for [to,ti] . (5.2) 

We consider a solution ^ of the Dirac equation with interaction (11. ip . Ac- 
cording to (j5.2p . ^'(t,x) is, for t < to, a solution of the free Dirac equation. 
We uniquely extend this free solution to the whole Minkowski space and 
denote it by ^in, i.e. 

(i^-m)§in = with §in(t,x) = §(t,x) for t < to- 

Similarly, ^'(t, x) is also for t > ti a solution of the free Dirac equation; we 
denote its extension by ^out) 

(z^-m)§out = with §out(i,^) = ^{t,x) for t > ti. 

The wave functions \I'in and ^'out are called the incoming and outgoing scat- 
tering states. The S'-matrix S maps the incoming scattering states into the 
corresponding outgoing states, i.e. 

§out = S §in for every § with (i^ S - m) § = 

As a consequence of the Dirac current conservation, 5" is a unitary operator 
(with respect to the scalar product ()1.2p ). Using the scattering states, one 
can build up asymptotic Dirac seas for t < to and t > ti. Namely, for an 
observer in the past t < to, the bosonic potentials are zero. Thus it is natural 
for him to describe the vacuum with the free Dirac sea (jl.3p . If this Dirac 
sea is extended to the whole Minkowski space with external field, one gets 
the object 

where the wave functions are the solutions of the perturbed Dirac equa- 
tion whose incoming scattering states are the plane wave solutions , 



{i^ + B- m) = with 
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Accordingly, an observer in the future t > to describes the vacuum with the 
Dirac sea 



a=l,2 



where 



{^i^ + B-m)^l = with (^y)„,t = *, 



ka " ^ ka^"""^ ka 

tion series 



The states ^'M and have a more exphcit form in terms of the perturba- 



oo oo 

= E(-«'^r^^a and *L = E(--"^r^.^a 

n=0 n=0 

as is immediately verified with (j5.2p using that the support of the advanced 
and retarded Green's functions is the future and past light cone, respectively. 
The asymptotics of the Dirac seas is completely described by the 5-matrix; 
namely 

^in = ^ = ^out (5.3) 

Po\= SP:^,S-' (5.4) 

-Pin = S ^hi ■ (5-5) 

The physical scattering process is conveniently described with the two Dirac 
seas of the observers in the past and in the future: If the physical system 
is described by P^, for example, the observer in the past is in the vacuum. 
According to ([5|) , does in general not coincide with the Dirac sea . This 
means that for the observer in the future, both positive frequency states are 
occupied and negative frequency states are unoccupied, so that for him the 
system contains both particles and anti-particles. This explains the physical 
effect of pair creation. Other scattering processes are described similarly. 

The causal perturbation expansion yields a unique object P describing 
the Dirac sea in the scattering process. P coincides neither with nor 
with P^; since its construction involves both the advanced and retarded 
Green's functions, it can be considered as being an "interpolation" between 
P^ and P^. At first sight, it might seem strange that the Dirac sea is now 
in both asymptotic regions t < to and t > ti described by the same object. 
Namely, it was essential for our discussion of pair creation that the Dirac 
seas of the past and future observers were different. It might seem that 
by redefining the Dirac sea, we no longer have pair creation. Clearly, this 
is not the case; all physical effects occur in the same way regardless if one 
works with the asymptotic Dirac seas P^, or with P. This is because 
the S'-matrix, which completely describes the physical scattering process, 
does not depend on the definition of the Dirac sea. Thus the choice of the 
definition of the Dirac sea in the asymptotic regions is merely a matter of 
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convenience. This may require some explanation: Suppose that we describe 
the Dirac sea with P. Then the asymptotic Dirac seas Pin and Pout consist 
of both positive and negative frequency states. As a consequence, they are 
not stable; they tend to decay into the Dirac sea P of all negative-energy 
states (this is clear physically from the fact that P has lower energy than 
Pin and Pout)- Taking this into account, one gets a consistent description 
of the physical observations. A further complication with P is that the 
current and energy distributions in the asymptotic regions are in general 
not homogeneous. For these reasons, it is highly inconvenient to describe 
the scattering process only with P; it is much easier to work with P^ and 
P^. But apart from these purely practical considerations, there is no reason 
against the description of the Dirac sea with P. The great advantage of 
the causal perturbation expansion is that it gives a unique definition of the 
Dirac sea, even in the region with interaction < t < ti. The Dirac sea 
is not defined with reference to an observer, but becomes a global object of 
space-time. 

Our definition of the Dirac sea is the starting point for the more technical 
analysis in [6], where all operator products are estimated and computed 
explicitly in an expansion around the light cone. In order to further clarify 
the definition of the Dirac sea, we now qualitatively anticipate some results 
of i. 

First of all, we explain what "causality" of the perturbation expansion 
for the Dirac sea precisely means: The expansion (j3.9p for km{x,y) is causal 
in the strict sense that the perturbation operator B{z) only enters for z in 
the "diamond" z G (L^ H Ly) U (L^ fl Ly). Since Pm{x,y) does not vanish 
for space-like y — x, its perturbation expansion, and consequently also the 
expansion of the Dirac sea, cannot be causal in this strict sense. As is 
shown in [6], the distribution P(x,y) has singularities on the light cone (i.e. 
for [y — x)^ = 0). It turns out that these singularities can be completely 
described in terms of B{z) and its partial derivatives along the convex line 
z G xy. Our perturbation expansion is causal in this weaker sense. It is even 
uniquely characterized by this "causality" of the singularities on the light 
cone. 

Both the operator products and the perturbation series were only treated 
as formal expressions throughout this paper. We outline in which sense these 
expressions make mathematical sense: It is shown in [6] that all operator 
products are well-defined distributions if reasonable regularity conditions on 
B are assumed. The convergence of the perturbation expansion is a more 
difficult problem. For chiral and scalar/pseudoscalar potentials, convergence 
is shown in [6j for the formulas of the light-cone expansion by explicit calcu- 
lation. For a gravitational field, the situation is more complicated, because 
the contributions to P(x, y) of higher order in B become more and more 
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singular on the light cone. With a Taylor expansion of the (^-distribution 

5{x + a) = 5[x) + a5'{x) + —5"{x) + ••• 

these contributions can be understood as describing a "deformation" of the 
light cone (corresponding to the diffeomorphism invariance of General Rel- 
ativity), but the convergence has not yet been established rigorously. 

We finally remark that the fermionic projector y) of section U] is 
considered in [3| as the basic physical object. In this context, the above 
construction gives a unique characterization of P by a perturbation B of the 
Dirac operator. This makes it possible to get a connection to the description 
of the interaction with classical potentials. It turns out that this "classical 
limit" is completely determined by the singularities of P{x,y) on the light 
cone. The "causality" of our perturbation expansion is then directly related 
to the locality and causality of the classical field equations. 



References 



[1] J. Bjorken, S. Drell, "Relativistic Quantum Mechanics," Mc Graw-Hill, 
1964 

[2] B. Thaller, "The Dirac Equation," Springer Verlag, 1984 

[3] F. Finster, "Local U{2, 2) Symmetry in Relativistic Quantum Mechan- 
ics," preprint hep-th/9703083^ to appear in J. of Math. Phys. 

[4] F. Finster, "Derivation of Field Equations from the Principle of the 
Fermionic Projector," preprint gr-qc/9606040| 

[5] A. Ashtekar and A. Magnon, "A Geometric Approach to External Po- 
tential Problems in Quantum Field Theory," Gen. Rel. and Grav. 12, 
205-223, 1980 

[6] F. Finster, "Light-Cone Expansion of the Dirac Sea in the Presence of 
Chiral and Scalar Potentials," preprint hep-th/9809019 



